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Abstract. We briefly review the strategy to perform non-perturbative heavy quark effective 
tlieory computations and we specialize to the case of the b quark mass which has recently been 
computed including the 1/m term. 



1. Introduction 

Lattice QCD allows a first-principle study of non-perturbative properties of QCD (e.g. 
computation of the hadron spectrum, decay constants and other matrix elements). However 
many systematics have to be controlled in order to reliably compute quantities of 
phenomenological interests: 

(i) numerical computations can be performed only at finite non-zero lattice spacing. Results 
thus have to be extrapolated to the continuum limit; 

(ii) the presence of an ultra-violet cut-off (the lattice spacing a) and of an infra-red cut-off (the 
volume V) constrains quark masses, and extrapolations are needed to reach the chiral and 
the heavy quark regimes; 

(iii) dynamical light quark effects are numerically expensive to simulate and are thus often 
neglected in the so-called quenched approximation. This approximation turned out 
empirically to be a good approximation in many cases and is useful in order to develop 
new methods and to pin down systematics. All the results presented in this contribution 
are obtained in this approximation. 

We now discuss the second point above in more detail. Pions have a Compton wavelength 
which is too large compared to lattices that can be presently simulated and in order to keep 
finite volume effects under control one has to simulate light quarks (u, d) which are heavier than 
the physical ones and then extrapolate to the physical point by matching lattice results with 
the chiral effective theory. 

On the contrary, the B meson is too heavy and its Compton wavelength is too small compared 
to the lattice spacings that can be simulated at present (for large enough volumes). In order to 
keep lattice artifacts under control one is forced to simulate quark masses in the region of the 
charm quark and then extrapolate to the b-quark mass. Another interesting possibility is to use 
heavy-quark effective theory (HQET) to describe the b quark. HQET can be formulated on the 



^ based on a work in collaboration with M. Delia Morte, N. Garron, R. Sommer and part of the ALPHA Coll. 



lattice as shown in p^. However, in order to be able to perform the continuum limit, it has to be 
renormalised non-perturbatively. A strategy has been devised and applied to the determination 
of the b-quark mass in the static limit of HQET [2] and presented in a more general framework 
in [3]. We will briefly review this strategy and then describe its application to the case of the 
computation of the b-quark mass including 1 /mb corrections ^ . 



2. Non-perturbative HQET 

The HQET action on a lattice (up to terms of order 1/m?) reads 

-^HQET = a^y]{V'h-^oV'h +t^spin'0h(-f^ • B)^h+WkmV'h(-5D^)V'h} +0{l/w?) (1) 

^ ^ — ;? ' ^ — ' ' 

where the first term is 0(1} [1] while the second and third terms are 0(l/m) {iVs-pm and Wkin 
are formally of order 1/m) q The heavy quark fields are subject to the constraints -P+V'h = V'h; 
4>hP+ — V'h with P+ = ^(1 + 7o). In the following we will consider two different discretisations 
of the heavy quark action where the naive parallel transporter in the covariant derivative Dq is 
replaced by HYPl and HYP2 gauge links (see [6] for details). 

Also the composite operators have an 1/m expansion in the effective theory. For example, 
the time component of the axial current, including 1/m terms, reads: 



yl^Q^^(x) = Z^Q'^^V^iToTs^h + c^ViT.^D.V'h +0(l/m2) (2) 

where c^^^ = 0(l/m). 

In the path integral one keeps the 0(1) term of the HQET action in the weight and 
expands the rest in powers of 1/m considering the higher order terms of the HQET action 
{C'^'^^x) = 0(l/m'^)) as operator insertions: 

{0)=Z-^ J D0e-^''^ht-a4 2^^hWSot/'hW Q |l-a^^/:(^)(x) + ...| (3) 

In this way the effective theory is renormalisable order by order in 1/m (this statement is 
equivalent to assert the existence of the continuum limit due to universality) and results in an 
asymptotic expansion in 1/m. It is important to notice that these properties are not automatic 
for an effective field theory: e.g. ChPT shares these properties while NRQCD does not. 

It is important to mention that for HQET on the lattice, as for any other theory with a 
dimensionful (hard) cut-off, the possibility of mixing with operators of lower dimension implies 
that renormalisation has to be carried out non-perturbatively, in order to be able to perform 
the continuum limit. 



3. Matching between QCD and HQET 

The bare couplings of HQET (mbare; ^kim "^spin, c^^^"^ , Z^^^"^ , . . .) are unknown parameters and 
have to be determined by matching the effective theory with QCD, thus allowing a transfer of 
predictivity from QCD to HQET. A way of determining the bare couplings of HQET in a non- 
perturbative way (thus assuring the existence of the continuum limit) is to perform the matching 
non-perturbatively. In practice this consists in properly choosing A'^hqet (the number of HQET 
bare couplings to be determined at a certain order in 1/m) QCD observables and to equate 
them to their corresponding expansion in HQET: 

^QCD^^HQET^ k=l^2,...,NnQET (4) 



^ for the unexplained notation we refer to [4] 



The bare couplings of HQET can consequently be extracted by solving the resulting system 
of equations. It is important to notice that this procedure requires to be able to simulate 
the b-quark around its physical mass. How is this possible? The trick consists in performing 
the matching in a small physical volume and at small lattice spacing (using the Schrodinger 
Functional) in such a way that mbO <C 1. At the same time, in small volume the expansion 
in 1/m will have corrections in \/m\yL and thus L can not be chosen too small. We choose 
L = Li K, 0.4fm which allows both conditions mta -C 1 and l/{m\yL) <C 1 to be satisfied. 

However, physical observables need a large volume, such that the S-meson fits comfortably, 
say L ~ 4Li ~ 1.6 fm. The connection between these different volumes is achieved by a recursive 
finite size method based on the concept of step scaling functions (SSFs) [5]: 

$^Q'^^(2L) = F, ({<i>™(L), J = 1, . . . , iv}) (5) 

All these steps allow a fully non-perturbative formulation of HQET where the continuum limit 
can be taken in all steps [3J. 

4. Mb static (at order 

We present the strategy described above in the simple case of Mb in the static approximation 
(i.e. at order l/wP). We start from the definition of the finite volume i?-meson "mass" 
r = — 9o log[/A(a;o)]a;Q=L/2,T=L which is defined in term of the Schrodinger functional correlator 
/a(3;o) of the temporal component of the axial current (inserted at time xq) with a pseudoscalar 
source on one boundary (the temporal extension T of the Schrodinger functional is here chosen to 
be equal to the spatial size L). From F we construct the observable <^2{L,M)^'~^^ = LT{L, M). 
In the effective theory this observable has the expansion ^2{L, M)^^^"^ = L[T^^^^{L) + mbare]) 
where 772-bare is the overall energy shift between the effective theory and QCD. The matching 
condition ^2{L, M)^'~^^ = ^2{L, M)^^^"^ in infinite volume becomes ttt-b = -Estat + "^bare where 
-E'stat = lim^^oo F'^*''*(L). Multiplying by L2 = 2Li and eliminating mbare one obtains: 

me = L2^stat - LaF^t^H^i) + ^^"^^(i^i, Mb) 



L2E,t.t - L2r'^HLi) + :^cI.QCD(Li, Mb) 



Li 



^2i?stat - L2F^*^*(L2) +L2r''^\L2) - LsF^*^* (Ll) +^<I>^^° (Li , Mb) 



= L2[Estat-F^^"^(L2)] + a^iui) + 2Lir(Li,Mb) (6) 

a^OinHQET a^OinHQET a->0 for A/b-Li>l 

where ara{g'^{Li) is the SSF for the static effective mass F'^*'^*. The whole procedure is represented 
by the following diagram: 
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Figure 1. example of continuum limit extrapolation for the QCD observable <&^^^ (for three 
values of the mass around Mb) and for the static SSF fTm- 



After having taken the continuum limit in the various parts of the last line of Eq. [6] one 
has to solve it for Mb (the RGI b-quark mass). In the simulations we fix Li ~ 0.4 fm (from 
(f'{Li) ~ 3.48), we set the light quark mass to zero and we fix the RGI quark masses of the 
heavy quark around Mb. In infinite volume (Lqo = 4Li) we set the light quark mass to the 
strange quark mass because as physical input we use the mass of the meson. 

Two examples of continuum limit extrapolation are shown in Fig. [TJ The solution of 
Eq. [6] is represented graphically in Fig. [21 where the red square represent the value of 
L2m-Q — L2[-Bstat — r'^*^*(Li)] — o"jn(5^(ii)) (with m-Q set to the experimental value) while the 
green circles correspond to the values of 2$^^^(Li,M) for tree values of z = LiM around 
-LiMb. From the solution of Eq. [6]one obtains also the slope 

S = -^^f^^°(Li,M) = 0.61(5) (7) 

which will be needed in the computation of the 1/m corrections. The results for Mb in the static 
approximation are 

M^*''* = 6806 ± 79 MeV (8) 
and the error is dominated by the error on the renormalisation constant of the quark mass [4] . 

5. Mb at order 1/m 

At order 1/m, beyond the parameter mbare one has to determine at least (depending on the 
strategy) the bare couplings Wkin and ujspm which are of order 1/m. In our case we consider the 
spin averaged i?-meson mass and in this quantity Wspin cancels out. We rewrite Eq. [6]in the case 
1/m: 

OO volume = £;stat _|_ mbare + l^kinE^^^ 

Matching 1 $^^^°(L,M) = u;],inRf''{L) = ^f^^^ 

Matching 2 ri(L,M) = Tf^\L) + mbare + ^kinrf"(L) = 



HQET 



where now ^^"""^ = LTi and Fi is defined similarly to F but in terms of the Schrodinger 
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Figure 2. Graphical solution of Eq. [6]to determine M^^^. 

functional boundary to boundary correlators (this allows a determination of Mb at order 
1/m without need of determining c^^"^) while Tf^^{L) and r^'°(L) arise naturally from the 
expansion of Fi in the effective theory. is a more complicate quantity whose definition 

(together that of i?^™) can be found in [1]. The important thing here is that its expansion in 
the effective theory is proportional to Wkin, thus allowing to eliminate this parameter from the 
above system of equations. By eliminating also mbare one obtains the following equation 



istat 



rf^'^(L) +Ti{L,M) + 



1.Q^°(L,Af) . kin 



rk-(L)) 



(9) 



In this equation we set L = L2 = 2Li and we use suitable SSFs to relate $j(L2) with <I>j(Li): 

^i{2L,M) = af''{u)'^i{L,M) , $2(2^, M) = 2$2(^, M) + cji^('u) + a^'^(«)$i(L, Af) (10) 

The continuum SSFs are defined in terms of those at finite lattice spacing as ct(m) = 
lim^/z^^^Q S(u, a/L) where the definition of the S's (neglecting their 9 dependence, see [3] for 
further details) is: 

i?f"(2L) 



-ikin. 



u, a/L) 



Em(n,a/L) 



2L [Ff °(2L) - Tf^'iL)] 



2L 
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7713'^* + mg"-* +mj^^\ M^*'^* is the solution of 



The r.h.s. of Eq.[9]can be split in three parts me 

^gat^jy^stat^ _ ^^av -^j^jig ^]^g 1/m correction can be found by solving the equation mB(M^*'^* + 

^ [m^ (M^t^t ) + m'i^'' (M^*^t )] 



M^^"^ + M^^''^) = mf , which translates into M^^"^ + M^^''^ 



From Eqs. [U]and[TU]one obtains the formulae for m^'^\ m^'^^: 

m^^''\M) = ^4-(ni)cI>i(Li,M) ^(M) = i^^^:^M^q>,{L2, M) (14) 

Some example of continuum extrapolations at order 1/m are shown in Fig. [3l These are the 
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Figure 3. Continuum limit extrapolations for mg"^ (combining HYPl and HYP2 actions) and 
for m^^^ (only HYP2 shown). 



most difficult steps of the computation because in both S^'" and {E^^"^ — T^^{L2)) there are l/a? 
power divergences that have to cancel out and the extrapolation to the continuum limit is only 
linear in a (in the case {^E^^^^ — r|*'^*'(L2)) the divergence goes only as 1/a and it extrapolates - 
as is the case for (^^^^ and - quadratically in a). The 1/m correction to Mb are then 
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and in the MS scheme: 
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-rW°(L2))*i(L2,M) 
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= -25(13) MeV 
24(32) MeV 



(15) 



mb(mb) 
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mb + ml^^ 



4.347(48) GeV , 



m 



(1) 



-0.027(22) GeV , 



(16) 



which agrees with the value given by the Particle Data Group, despite the quenched 
approximation. 

An alternative computation has been performed by using observables defined in terms of 
/a(xo). In this case there is one additional combination of bare couplings of order 1/m that has 
to be determined and thus more step scaling functions are required. The final result is found to 
agree up to (small) 0(l/m-^) corrections. 
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